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A formalism for energy-dependent many-body perturbation theory (MBPT), previously indicated 
in our recent review articles (Lindgren et ai, Phys.Rep. 389,161(2004), Can.J.Phys. 83,183(2005)), 
is developed in more detail. The formalism allows for a mixture of energy-dependent (retarded) 
and energy-independent (instantaneous) interactions and hence for a merger of QED and standard 
' (relativistic) MBPT. This combination is particularly important for light elements, such as light 

, heliumlike ions, where electron correlation is pronounced. It can also be quite significant in the 

• medium-heavy mass range, as recently discussed by Fritzsche et al. (J.Phys. B38, 8707(2005)), 

' with the consequence that the effects might be significant also in analyzing the data of experiments 

with highly charged ions. A numerical procedure for treating the combined effect is described, and 
some preliminary numerical results are given for heliumlike ions. This represent the first numerical 
evaluation of effects beyond two-photon exchange involving a retarded interaction. It is found that 
for heliumlike neon the effect of one retarded photon (with Coulomb interactions of all orders) 
■ represents about 99% of the non-radiative effects beyond energy-independent MBPT. 
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' What is commonly known as Many-Body Perturbation Theory (MBPT) is a class of pcrturbativc schemes for 
] bound atomic, molecular or nuclear states with a time- or energy-independent perturbation, based upon the Rayleigh- 
"Y^- Schrodinger perturbative scheme, such as the Brueckner-Goldstone linked-diagram expansion or variants thereof [j]. 
Also certain iterative or "all-order" approaches, like the Coupled-Cluster Approach (CCA) or "Exponential Ansatz", 
can be referred to this category, although they are not strictly perturbative in nature. All these approaches can in 
' principle treat the electron correlation to arbitrary order. In addition, they have — as distinct from procedures based 
upon the Brillouin-Wigner (BW) perturbation expansion — the important property of being size extensive in the sense 
that the energy scales linearly with the size of the system. Furthermore, by using an extended or multi-reference 
model space such schemes can also successfully handle the quasidegenerate problem with closely spaced energy levels 
that are strongly mixed by the perturbation. 



For time- or energy dependent perturbations, like those of quantum-electrodynamics (QED), the situation is quite 
different and much less developed. There is presently no numerical scheme available that can treat energy-dependent 
perturbations together with electron correlation to arbitrary order, and also the treatment of quasidegeneracy forms 
a serious problem in connection with such interactions. In the numerical methods presently available for QED 
calculations the electron-electron interaction is treated by the exchange of fully covariant photons, which is quite a 
tedious — and usually unnecessarily tedious — process to handle the electron correlation. At most two-photon exchange 
can be treated in this way with computers available today, which is insufficient for light and medium-heavy elements, 
where the electron correlation is quite important. 



A major problem in extending the energy-dependent perturbation theory to include electron correlation is that 
most methods have a structure that is quite different from that of energy-independent perturbation theory, which 
makes it difficult to utilize the well-developed methods of the latter. Of the available methods only the Covariant 
Evolution Operator (CEO) method that we recently developed, has a structure that is akin to standard energy- 
independent MBPT [llllll. This opens the possibility of combining the two approaches as proposed in our recent 
review articles 0, and further developed in the present work. In that scheme the CEO method is combined with 
all-order MBPT methods of coupled-cluster type, so that the exchange of covariant photons can be mixed with an 
arbitrary number of instantaneous Coulomb interactions. In this way electron correlation can for the first time be 
treated to arbitrary order together with energy-dependent interactions of QED type. 
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The CEO method was originaUy developed in order to be able to treat the quasidegeneracy problem in QED 
calculations. The standard procedure for bound-state QED is the S'-matrix formulation 0, which has been successfully 
applied particularly to highly charged ions. For lighter elements — in addition to the electron-correlation problem — 
also the quasidegeneracy problem might be quite pronounced, and the S'-matrix formulation fails. One illustrative 
example is the fine structure of heliumlike ions, where in the relativistic formulation for instance the lowest triplet 
state ^Pi is a mixture of the basis states ls2pi/2 and ls2p3/2, which are very closely spaced in energy and strongly 
mixed for light elements. In order to be able to use the procedure of standard MBPT, where quasidegenerate states are 
included in an extended model space, off-diagonal elements of the effective Hamiltonian have to be evaluated, which 
is not possible in the standard S'-matrix formulation, due to the energy conservation of the scattering process 0- In 
the CEO method, based upon the evolution operator for finite time Q , the extended- model-space technique can be 
used, and a few years ago we applied this to evaluate QED contributions to the fine structure for heliumlike ions, 
including the quasidegenerate ^Pi state, down to Z = 9 Q. The same problem has very recently been treated by 
Shabaev et al. down to Z = 12, using the two-times Green's function technique Furthermore, the presently 

available methods for numerical QED calculations suffer from the shortcoming that they are not applicable to the 
lightest elements — below Z=10, say — due to convergence problems. The extension of the CEO method to include 
electron correlation to arbitrary order, presented here, is expected to remedy this problem. 

The most accurate calculations on ligh t heliumlike ions have been performed by the analytical or " unified" method 
of Drake and related techniques 0, 

Hi III III El- This is based upon expansion in powers of the fine-structure 
constant a of the Bethe-Salpeter equation and the Brillouin-Wigner (BW) perturbation series 0, 0| . It leads to 
extremely high accuracy for the lightest elements, but there is still a significant discrepancy between theory and 
experiment for the lowest triplet state of neutral helium. Being based upon the BW perturbation expansion, the 
procedure is not size-extensive and less suitable for larger systems. 

In addition to the light heliumlike ions and other light systems, the combination of QED and many-body effects 
can be of importance also for heavier systems, as recently discussed by Fritzsche et al. [T^ . They have analyzed 
this problem particularly with regard to possible heavy-ion experiments of the type that can be performed at the big 
storage rings, like that at GSI in Darmstadt. They then conclude that the accurate treatment of the interplay between 
QED and many-body effects constitutes one of the most challenging problems in connection with highly-charged-ion 
experiments. 

The present paper will be organized as follows. In the next section we briefiy review the standard perturbation 
theory for time-independent and time-dependent perturbations, and next we summarize the properties of our recently 
introduced covariant evolution-operator method. The following main section deals with the derivation of Bloch 
equations for combined retarded and instantaneous interactions, which constitute our working equations for treating 
this problem. Finally, we describe briefly our numerical procedure and give some preliminary numerical results, 
including the first numerical result of effects beyond two-photon exchange with a retarded interaction. The numerical 
procedure together with more complete numerical results will be published separately |l8| . 



II. STANDARD MANY-BODY PERTURBATION THEORY 



A. Time-independent perturbation theory 



In the multi-reference form of MBPT we consider a number of target states that are eigenstates of the Hamiltonian 
of the system 

H \<f") = \^") ; (a = l,2,---d). (1) 

The Hamiltonian is partitioned into a model Hamiltonian, Hq, and a time-independent perturbation. H' , 

H = Ho + H'. (2) 

For an iV-electron system the model Hamiltonian is assumed to be composed of single-electron Schrodinger or Dirac 
Hamiltonians 

N 

^o-^/ioW- (3) 

i 

For each target state there is a model state, confined to a subspace, the model space, with the projection operator P. 
Tn the intermediate nnrma.lizatinn we ime. here the model states are the nrniertinn of the cnrresnnndincr tarcret state.s 
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on the model space 

l*o) = ^^1*")- (4) 

A wave operator can be defined for the inverse transformation 

\^n = n\^^); (a = l,2,-..d). (5) 

An effective Hamiltonian can be defined, Hcs = PHflP, that operates in the model space and for which the eigen- 
vectors are the model states and the eigenvalues the corresponding exact energies 

i/eff|^'?)=i;"|*^). (6) 

The corresponding effective interaction is defined 

Kff = i?cff - PHoP = PH'nP (7) 
The wave operator satisfies the generalized Bloch equation 0, |^ 



n,Ha]P = {H'n -QVeff)P (8) 



This leads to the Rayleigh-Schrodinger perturbative expansion for a general multi-refcrence (quasidegenerate) model 
space, and it can also be used to generate the corresponding linked-diagram expansion of Brueckner-Goldstone type 

[n,Ho]p = {H'n-nv,^)^^^^^^^p (9) 

Here, only so-called linked terms or diagrams survive on the right-hand side. Using the exponential Ansatz 

n = {c^}, (10) 

where the curly brackets represent normal-ordering, leads to the Coupled-Cluster expansion for the same model 
space 

[5,i/o]P=(i/'f^-f^Kff)_,P. (11) 

Here, all terms on the right-hand side are "connected". (For the distinction between "linked" and "connected", see, 
for instance, ref. 0. 



B. Time-dependent perturbation theory 



In this paper we shall only be concerned with stationary states, but we need for our purpose to use the formalism 
for time-dependent perturbation, which we shall briefly review. 

The time-dependent state vector satisfies the time-dependent Schrodinger equation (using relativistic units h = c = 
m eo = 1) 

il\^it))^Hit)\xit)}, (12) 

where 

H{t)^Ho + H'{t). (13) 

Hq is the time- independent model Hamiltonian and H'{t) is a perturbation that might be time-dependent. In the 
interaction picture (IP), where an operator is related to that in the Schrodinger picture (SP) by 

Oi(t) = e'^"*Ose-'^''* (14) 

the Schrodinger equation becomes 

i^|Xi(t))=i/i'(<)|xi(i)). (15) 

The time-evolution operator is defined by 
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and satisfies the equation 



with the solution |22 



i^Ui{t,to)^Hi{t)Uiit,to) (17) 
ot 



(18) 



Uiit,to)^l + J2^-r d^^n--- / d''xiTu[n{{xM{xn-i)---H[ixi)]. 

Here, Td is the Dyson time-ordering operator, and Ti.'i{x) is the perturbation density, defined by 

H[{t) = J d^xn{{t,x). (19) 

An adiabatic damping is added to the perturbation 

H[{t) ^ = H{ e-^l*! ; C/i(t, to) "> f/7(i, ^o), (20) 

where 7 is a smaU, positive number. This imphes that as t — s- ±00 the eigenfunctions of H tend to eigenfunctions of 
Ho. 

For stationary states we assume that the perturbation H' is time-independent in the Schrodinger picture — apart 
from the adiabatic-damping factor. Then according to the Gell-Mann-Low (GML) theorem ^2^1 the wave function at 
time t = is for a single-reference model space given by 

l^)Hx(0)) = linr ,;^(°'-)|^°) , (21) 

where \'^o) is the time-independent model state. This can be generalized to a general multi-reference model space ^ 
Eq. 110] 

|vl.-)^lim,^"f^--("'— («=l,2,...d), (22) 
' ' 7-0 ($"|i7^(0,-oo)|$") ' ^ ' ' ^' V ; 

where is a normalization factor and the vector |$") is defined 

|$") = lim lim \x°'{t)). (23) 

7^0 t—* — oc 



Then the wave function at time t ~ satisfies the time-independent Schrodinger-like equation 



(24) 



where H' is the time-independent perturbation in the Schrodinger picture. 



The evolution operator normally contains singularities or quasi-singularities as 7 ^ 0, when an intermediate state is 
degenerate or closely degenerate with the initial state. In the GML formulas these (quasi-)singularities are eliminated 
by the denominator so that the ratio is always regular, which is one formulation of the linked-diagram theorem [2^ . 



III. COVARIANT EVOLUTION-OPERATOR APPROACH 



(The reader is referred to refs |j, |3 for more details concerning the basic covariant evolution-operator formalism.) 



A. Interaction with the electro-magnetic field 



We consider now the interaction between electrons and the quantized electro-magnetic field represented by the 
perturbation density 
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Here, '4>i, V'l are the electron field operators, the Dirac alpha operators and is the radiation field 

oc (al{k) e*'^" + a.(fc) e"'''^) , (26) 

applying standard summation convention, aj(fe), ar{k) are the photon creation/absorption operators, k is the wave 
vector, K, is the four-vector momentum n = (w, — A;) (fc = |A;|), and represent the polarization vectors. The 
Hamiltonian of the radiation field is represented by 



-ffrad = UJ a\.{k)ar{k) = k al{k)ar(k). 



(27) 



t^t' 



FIG. 1: Graphical representation of the covariant-evolution operator for single-photon exchange in the form 13m (left) and in 
the form I31II with to — * —oo- 



With the interaction density H25|l the evolution operator (|18f) for single-photon exchange becomes 

integrated over all space coordinates and time coordinates as indicated. ipl^ are the positive-energy part of the 

electron-field operators, and the interaction kernel 

I 1 

iI{xi,X2) = {-ea^Afj,)i{-ea'^A^)2 = e'^a'^a2\DF^j.v{xi - X2) (29) 

is given by the product of two perturbations H25|) with contraction of the radiation-field operators (indicated by the 
hook). D^^v{xi — X2) is the Feynman photon propagator. 



The evolution operator above is non-covariant but can be made covariant by inserting zeroth-order Green's functions 
on the in- and outgoing states 0, 0| 

Uclit'^to) = '^''x',A^x'^^l{x\Wi{x'2) 11 d^x^d^X2Go{x[,x'^;xi,X2) 

X J J d^a;iod^a;2oi^(a:^i,2:2)Go(a;i,a;2;a:io,a;2o)'/'i(a;2o)'^'i(a;io)- (30) 

Here, the time integrations over ti and t2 are performed over all times, and positive- as well as negative-energy states 
are allowed as incoming and outgoing states. The initial and final times arc the same for the two electrons, i.e., 

^10 ^ ^20 ^ ^0 and t-^ — t^ = t . 

We shall in the following assume that the initial time is ~ ~oo, and then due to the adiabatic damping l|2()(l we 
can leave out the rightmost Green's function Q 

t^Cov(^''~oo) = ^i^l{x[)-4'l{x2) Gq{x[, x'^; xi, X2) il{xi, X2) il^i{x2)4>i{xi). (31) 

Here, we have left out the integrations, and in the following we shall also leave out the subscript cov as well as the 
initial time. 



B. Wave operator and effective interaction 

The evolution operator is generally singular and can be expressed 

n(t) P = P4- U(t.) P . PTI(O) P. (^2) 



FIG. 2: Graphical representation of single-photon potential 1411 1. 



where U{t) is always regular and known as the reduced evolution operator — all singularities are collected in the 
last factor PU{0)P. The heavy dot indicates here that the two factors evolve in time independently from different 
model-space states. For the time t = this becomes 



U{0)P = [l + QU{0)]P- PU{0)P, 
where Q is the projection operator for the "complementary space" (outside the model space) 

Inserting the expression H33|) into the GML formula H22|) . yields 

1^") . [1 + g^(o)]P hm ^ [i + Qm^i^n- 



(33) 



(34) 



But = vJ/q is the model state and hence the expression in the square brackets represents the wave operator 



n = i + Qu{o) 



(35) 



The Q operator is for a two-electron system given by 

Q = l-P= \rs){rsl 



(36) 



where the ket vectors \rs) represent straight (non-antisymmetrized) products of single particle states, summed over 
all states outside the model space. The single-particle states are generated by the single-particle Hamiltonians (0) 



(37) 



which also include the nuclear field (Furry picture). In the two-electron system we shall study here, the occupied 
electron states are treated as open-shell or valence states, which implies that there are no core or hole states — apart 
from the negative-energy states. 



Performing the integrations for the single-photon exchange (|31f) . yields 

{rs\U^^\t)\ab) = (rs|e-"(^--^'-^»)rQ(£;,6)Vi(i?„fc)|a5) 
or, generally, operating on a model-space state of energy 

U^'^\t)P = e-'*(^-^°'rQ(f )Vi(f )F. 

Here, 

Q 



Eab — Sa+ £b 



rQ(^) = 



is the resolvent, and 

{rs\Vi{£)\tu) ~ Irs I dkf{xi,X2,k) 



£ - Er - £u - {k - i7)r £ - Es - et - {k - i'y)s 



tu 



(38) 
(39) 

(40) 
(41) 



is the matrix element of the potential, considering both time orderings (see Fig. I^J. The subscript represents the 
sign of Er- Using the evolution operator (|39|) and the relation H35|) (with f/'^' = J/'^-*), this gives the wave operator 
for single-photon exchange 



n^^^P = rr.(£)V-,(F.)P. 
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The effective interaction ^ can generally be expressed 0| 



P 



(43) 



The reduced evolution operator has the same time dependence in all orders as in first order H,S9|I . which implies that 
the time derivation eliminates the resolvent. In first order this yields 



V^^{£)^PV,{£)P. 



(44) 



The function /(a^i, k) in the potential (|41|) depends on the gauge used and is in the Fcynman gauge given by 

r r ,N sin(fcri2) 

fF(xi,X2,k) = -— ^ (1 - cci • a2) , 

47r"^ ri2 

where ri2 is the interelectronic distance. 

In the Coulomb gauge, which is natural to use in many-body calculations, the potential can be separated into an 
instantaneous and a retarded part, 



V{£)^Vi + VR,t{£), 

where only the latter is energy dependent. The instantaneous part is the Coulomb interaction 
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and the retarded part is given by the expression ()41|l with 



fc{xi,X2,k) = — 



sin(fcri2) , ^7 ^^^^^2) 

Otl • OL-2. V \OL\ • Vl) (0:2 • V2 



r\2 



ri2 



(45) 



(46) 



(47) 



where the nabla operators do not operate beyond the square bracket. Here, the first term represent the Gaunt part 
and the second term the scalar-retardation part, which together form the Breit interaction. In the following wc shall 
assume that the Coulomb gauge is used. 

The relations given here, particularly the framed equations (|35|l and H43|l . demonstrate the close analogy between 
the covariant evolution-operator approach and standard MBPT, which opens up the possibility for a merger of the 
two procedures. 



IV. BLOCK EQUATION FOR INSTANTANEOUS AND RETARDED INTERACTIONS 

A. Retarded interactions 

From the definition H32|) . the following relation can be derived, using standard algebra |^, 

U{t)P = U{t)P + U{t) {PUP - -PUP), 



where 



U{t)P = U{t)P - M{t) 



(48) 



(49) 



is the evolution operator without intermediate model-space states and M is the model-space contribution (MSC). The 
JJ operator in Eq. H48|) has the time argument t = 0. 

As mentioned, in the product U {t)PU (0)P the time-dependent part U (t) evolves from the energy of the state to the 
far right (f ), while in the dot product U{t) • PU{0)P the operator U{t) evolves from the energy of the intermediate 
state {£')■ Therefore, the second term, which represents the MSC part, becomes 



M = (TKF.) - TKF.'W PTIP. 
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This can be expressed 



PUP = 



off 



where is the analogue of the effective interaction (|43|l without MSG 

d 



P 



P. 



This yields 



U{£)-U{£') dU 



using the difference ratio defined in Appendix lEl This yields 



SU 



U{t)P=U{t)P+—VaS. 

o£ 



With the wave-operator relation (|35|l we then have 



(51) 



(52) 



(53) 



(54) 



where the last term represents the MSG (including "folded" diagram) and ft is the wave operator without MSG. 

The relation (|48|l leads to the expansion 

U{t)P = U{t)P + U{t) {PUP - -PUP) + U{t) {PUP - -PUP) {PUP - -PUP) + • ■ • (55) 

The exchange of a single photon corresponds to the second-order evolution operator and leads to the result given 
above H42|) . The two-photon exchange corresponds to the next even order of the evolution operator 

[7(4) {t)P = f7('') {t)P + [7(2) [t] {PU'-^^P - ■PU'-^^P) . (56) 

Since there is no MSG in lowest order, we have [/^^^ = [7'^) = U^^\ and the corresponding wave operator becomes 



5£ 

where fl^^^P ~ rgVirgViP. In the case of degeneracy this goes over into 



(57) 



(58) 



In third order we have App. D] 



(1) ^ -517(1) 

5£ ' 



which leads to 



d£ 



S£ 



d£^ 



(l)^2 



(59) 



(60) 



where the second-order difference ratio is defined in Appendix IbI The last term is associated with double model-space 
contributions ("double fold"). This leads to conjecture for the all-order expansion (c.f. Eq.116]) 



71=1 



(61) 



which we shall now verify. 
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In order to prove the relation above, we start by using this relation to form the difference ratio 



5£ 5E ' S£-2 '"^'SS 



n=l 



1 



SV.ff 



5£ 



(62) 



(It should be noted that the different V^s operators are in general associated with different energies, as explained in 
Appendix El Eq- IBSI ) Next, we take the time derivative of the relation (jHSJ, using the relation 143|) . which yields 



58 



This gives with the relation H62|l 



5£ 



(63) 



(64) 



and with the relation (|53|l we retrieve the relation H61|l . which is then proven. 

In order to obtain more Bloch-like relations, we first introduce the reaction operator, which is the effective interaction 
l)43|l ■ apart from the projection operators. 



This gives Vcft = PVrP, and the reaction operator is equal to the wave operator, apart from the resolvent, 

QU{0) = Qn = TgYn. 

We also introduce the operator 



where U is the evolution operator without MSG. This leads in analogy with the relation (|66|) to 
Using the rule for differentiating a product, developed in Appendix IbI CEa. IB4(I . and the relation 



j(m-l)rQ 



we can express the difference ratios of as 



The last term can be expressed 



E 

m=0 



6£ 



(65) 



(66) 



(67) 



(68) 



(69) 



(70) 



m=0 



Q 



(71) 



Inserted in the expansion (|61|1. this yields 



(72) 



This equation will later be used to derive the Bloch equations for energy-dependent interactions. 
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FIG. 3; Graphic representation of pair functions pab generated by the pair equation (I75II . 



B. Instantaneous interactions 



The instantaneous Coulomb interaction can be treated essentially as in standard (relativistic) many-body theory, 
and we start by recalling the treatment of the correlation effect for a two-electron system. The wave operator in the 
coupled-cluster formalism (fTU)) can then be expressed 



where S2 is the two-body cluster operator. Operating on a model-space state, yields a pair function 

\pab)^^\ab) = \ab) + s2\rs). 
Inserting the pair function into the Bloch equation , yields the corresponding pair equation 

{Sa+Sb - /jo(l) - ^o(2))|Pa6> = \r s) {r s\Vi2\ p ab) " IPcrf) (cd| Kff |a6) . 



(73) 
(74) 
(75) 



The last term is the folded term and is the result of the reduction of singularities which appear when the intermediate 
states lie in the model space, which we have referred to above as model-space contribution (MSG). When the equation 
is solved iteratively, the Coulomb interactions are generated to all orders, as illustrated in Fig. Bl This is the type of 
pair functions we have been using in our many-body calculations for several decades [23, |2^ l29l l30l l3ll | . 

We denote the wave operator with only Coulomb interactions by fJi and the part with no folded diagrams by fij. 
Then we have 



1+TqVi2+TqVi2TqVi2 



(76) 



and using the relations H72|) and this leads to the standard Bloch equation © 



(77) 



C. Combined instantaneous and retarded interactions 



We shall now find Bloch equations for the combined retarded and instantaneous interactions. In the Coulomb gauge 
the function f(xi,X2,k), involved in the exchange of a retarded photon, is given by the expression (|47|) . which can 
be separated into products of single-electron operators, as shown in Appendix 1X1 



k 

fc{xi,X2, fc) = ^ E [ - (2/ + l)V^ikr,) . V^{kr2) 



1 



1=0 



21 + 1 



Vinikri) ■ Vi^ikr2) 



(78) 



The two terms represent the Gaunt and scalar-retardation parts, respectively, and we shall treat each of them as 
the re.sult of tiiin nertiirha.tinn.'!. na.melv Vlikri) and Vj.ikv'i) in the c.a.'ie nf the Gaunt interaction an and 
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Vgp{(fcr2) for the scalar retardation — of course, with the appropriate factors and with summation over the angular 
momentum of the photon, I, and integration over space and the linear photon momentum, k. 

The photon can also be absorbed by the same electron, leading to self-energy and vertex-correction contributions, 
as we shall briefly indicate below. 

The perturbations above are time-independent in the Schrodinger picture, and we can then apply the Gell-Mann- 
Low theorem (|22ll . which leads to the Schrodinger- like equation (|24|) 



(79) 



where H' is given by the instantaneous Coulomb interaction Vi = V12 and the two retarded components VQ(fcr) and 
Vgj^(fcr). The wave function lies here in an extended Fock space with variable number of uncontracted, virtual photons, 
which we indicate by using the bold- face symbol The bold-face symbol Hq represents the model Hamiltonian 
including the radiation field H27|l . We also introduce a Fock-space wave operator in analogy with the standard wave 
operator © 



(80) 



The resolvent (|40|) is now generalized to 



Q 



S Ho 



(81) 



when operating on a model-space state of energy £. Q = 1 — P is here the projection operator for the complementary 
Fock space, for a two-electron system given by 



Q \rs){rs\ + \ij,k){ij,k\ 



(82) 



The first term represents the part of the operator in the restricted space with no photons (c.f. Eq. with \rs) being 
a state outside the model space. The second term represents the part with one photon, with \ij) being an arbitrary 
state, etc. 



b a 



b a 



V 
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+ 
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FIG. 4: Graphical representation the pair equations 19111 . solving the Bloch equation equations 18811 and 19Ut . 



In order to treat the case where the instantaneous interactions cross a retarded photon, we have to apply the former 
between the two perturbations of the retarded interaction. We denote the wave operator with an "uncontracted" 
retarded photon and an arbitrary number of instantaneous interactions before and after the retarded photon is 
created by ^Q{k) and f2gj^(fc), respectively, for the two components of the retarded interaction. The components of 
these operators with no model-space contributions are in analogy with previous cases denoted by ^^(fc) and rjgj^(A:), 
respectively. We then have for the Gaunt interaction (and similarly in the scalar-retardation case) 

n'Q{k)p = (1 + TqV,2 + TQVurQVi2 + ■ ■ ■)TQVl{k){i + TqVu + rQVuTQV,2 + ---)p. (83) 

The rightmost bracket represents the wave operator f2i (|76|l . which leads to 



fi'Gik)P = TQV^ik)niP + rQVi2n'cik)p. 



(84) 



12 





, s r , 




, s r , 




. s r , 




. s r , 






■ b a. 




. b a. 




+ 

. b a. 




■ b a. 





: H h folded 



< 



r , 




, s r , 








a, 




. b a. 



+ 



< 



H H folded 



FIG. 5: Graphical representation of a single retarded photon with crossing Coulomb interactions and without Coulomb inter- 
actions before and after the retarded interaction. 



Inserting this into the expression (|72|l . yields 



n=0 



I'' A/ ^" _L T 



(v;ff)"-rQr!'G(fc)\4ff, 



Since V12 as well as VQ(fc) are energy independent, we have 



n=0 ?i=0 

again using the expansion Ht)l|) and the rule ljB4p . Treating the second sum similarly, leads to 
and to the Bloch equation 



n'cik), Ho]p = v^{k)niP + Vi2^'Q{k)p - n'Q{k)v,s 



(85) 



(86) 



(87) 



In the next step the photon of the function ^\;^{k) is being absorbed by the other electron, followed by additional 
Coulomb iterations. We denote the corresponding wave operator by Via- Omitting for the time being the MSG 
associated with the Coulomb interactions, this leads to 



naP = (1 + TqVi2 + TqVv2TqVi2 



•TQV^{k)n'^{k)p, 



(89) 



integrating over k and summing over I according to the relation (|78|l . The MSC are as before obtained by inserting 
this relation into the formula (|72|) . which — using the same argument as before — yields 



^Ig..Ho\p = v^{k) n''^{k)p + Vi2nGP - naV.s 



(90) 



The framed equations above represent our main equations for dealing with the combined retarded and unretarded 
interactions. For a two-electron system they can be converted to the pair equations 

{ea + £b - ho{l) - h„{2) - k)\p'^^,,{k)) = \rs){rs\V^{k)\pab) + k-s) (rs|yi2 |p[,,,,(fc)) - \p'^^,^{k)){cd\V,^\ab) 

{ea + £b - /lo(l) - ha{2))\pG.ab) = \rs){rs\V^{k)\p'Q^,,{k)) + \rs){rs\Vi2\ PG,ab) — \pcd){cd\Vcs\ah) . (91) 

Note that the first of these equations is deduced from the relation (|88|l with the extended Hq, which leads to the 
momentum k on the left-hand side. These equations, which can be solved using standard technique |32l Is^ . are 
illustrated in Fig. 0] In Figs and El (upper line) we show more explicitly the diagrams involved in the process just 
described — in the former case without any Coulomb interactions before and after the retarded interaction and in the 
latter case with such interactions. When the photon is instead absorbed by the same electron as it is emitted from, 
we get the corresponding self-energy and vertex-correction effects — of course, after appropriate renormalization — as 
indicated in the bottom lines of the same figures. 
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FIG. 6: Graphical representation of a single retarded photon with crossing Coulomb interactions and with Coulomb interactions 
before and after the retarded interaction. The diagrams in upper line are obtained by solving the pair equations 1911 and those 
in lower line by analogous equation. 
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FIG. 7: Diagrams evaluated for the ls2s and ls2s^S states of heliumlike neon with the results given in Table n 



D. Derivation of the general Bloch equations 



The Bloch equations derived above are valid only in the case of no more than one uncontracted photon at each 
instance. In order to derive the more general Fock-space Bloch equation, we can proceed exactly as in the energy- 
independent case |23|, starting from the Schrodinger-like equation (|79|1 . We first project this equation on the model 
space 



•a -it, a 

' 



P{Ho + H')n^^ = Ho^o + PH'n^^ = E"^ 
using the fact that P and Hq commute, and then operate from the left with ft 

nHo^o + ^PH'n'^'i; = s"*". 

Subtracting the original SE H24|l . then yields the Fock-space Bloch equation 



n,Ha\P = H'nP-nV,s; H,s = PHoP + V,s; V,fi = PH'nP 



The solution can be expressed 



Qnp = TQ{H'n - nVcs)p 



(92) 
(93) 

(94) 
(95) 



where Tq is the generalized resolvent (|81fl . 

We can expand the Fock-space wave operator and resolvent into components acting in the subspace with no photons, 
with one photon etc. as 



n = f7 + r2+ + • • ■ 
rQ = rQ + r+... 



(96) 
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FIG. 8: Some irreducible diagrams, not included in the procedure described here. 
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FIG. 9: Examples of reducible diagrams, obtained by repeated use of the procedure described. 



The generalized Bloch equation 195() can then be separated into 



np = p + T, 



Q 



n+p = r 
n++p = 



Q 



Vi2n + VR,t^+ -nv. 



cff 



p 



cff 



p 



(97) 



The hook represents integration over the photon momentum k and summation over the angular momentum ^, according 
to the single-photon expression (|78|) . These equations are valid also in the case of multiple free photons at the same 
time. Considering at most a single free photon, we see that they lead to the Bloch equations (|88|) and H9U|) derived 
above. 



V. NUMERICAL PROCEDURE AND RESULTS 



The pair equations considered here can be solved numerically with essentially the same technique as dev elop ed by 
Salomonson and Oster for (relativistic) many-body calculations [3^ l3^ and used in our previous works 0, |3 ■ The 
radial integrations are performed with an exponential grid with 70-150 grid points and the k integration with 100-150 
points using Gaussian quadrature. For excited states poles appear in the k integration, which require special attention 
(for details, see ref. [3^). The numerical calculations are quite time consuming in the present case, since separate 
pair functions have to be evaluated for each value of the photon momentum. On the other hand, the procedure is 
particularly well suite for parallel computing, and we hope that the procedure can be speeded up considerably, when 
our routines are better optimized. 



TABLE I: Effects of one- and two-photon exchange for the excited ls2s and states of heliumlike neon (in fiH). 







ls2s^S' 


Is2s'-^S 


One-photon 


Gaunt 
Scalar ret. 


2465.44 
171,58 


171,50 
-171,58 


Two-photon 


Coul.-Gaunt 
Coul.-Scal.ret. 


-794.8 
22.5 


-51.8 
42.5 


One + two-photon 


Gaunt 
Seal. ret. 


1670.7 
194.2 


119.7 
-129.1 


One photon correlated 


Gaunt 
Seal. rot. 


1752.0 
183.9 


124.6 
-132.2 
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Here, we shall only give some illustrative examples of our numerical results — more results will be published sepa- 
rately [ig. In Tabic □ we show the effect of one retarded photon with one and infinitely many non-crossing Coulomb 
interactions for the ls2s ^5* and "^5 states of hcliumlikc neon. This corresponds to the diagrams shown in Fig. [7\ 
excluding virtual pairs (NVP). This represents the first numerical evaluation of effects beyond two-photon exchange, 
invol ving a retarded interaction. The two-photon effects have been compared with the corresponding S'-matrix re- 
sults |34| and are found to agree to 3-4 digits, which represents the numerical accuracy of the present calculations. 
The effect of correlation beyond second order is in this case found to be about five percent of the total contribution, 
which is one order of magnitude larger than the effect of the retarded (reducible and irreducible) two-photon inter- 
action |33- The effect of virtual pairs is of the same order as that of two retarded interactions. This indicates that 
for heliumlike neon the procedure described here with a single retarded photon together with Coulomb interactions 
represents about 99 % of the non-radiative effects not included in a standard many-body perturbation treatment with 
only instantaneous Coulomb interactions. For lighter elements the importance of a single retarded photon is even 
more pronounced. 



VI. SUMMARY AND OUTLOOK 



We have in previous articles described a new technique for QED calculations that we refer to as the covariant 
evolution-operator (CEO) method 0j S S I3- This method has the great advantage compared to other QED 
techniques that it has a structure very akin to that of standard many-body perturbation theory, which opens up the 
possibility for a merger of the two techniques. In two previous papers |j, [Sj we have indicated how such a merger 
can be accomplished, and this is developed further in the present paper, and some numerical results arc presented. 
Combined QED and correlation effects, which can be treated only in a very limited fashion by standard techniques, are 
of particular importance for light and medium-heavy elements. The CEO method also has the advantage compared 
to the standard S'-matrix technique that it can be applied to the case of quasidegeneracy (a property it shares with 
the two-times Green's- function technique of Shabaev et al. @). 

The procedure we have developed represents the exchange of a single retarded photon and an arbitrary number of 
instantaneous Coulomb interactions between the electrons, crossing and non-crossing. So far we have been working 
with positive-energy intermediate states — no-virtual-pair (NVP) approximation — but single and double virtual pairs 
can be included in the procedure. The procedure can also be used — with proper renormalizations — for radiative effects 
(self energy and vertex corrections) with a single retarded photon (see Figs IHl and [S)) . In principle, the procedure can 
be used also for irreducible multi-photon effects, where retarded interactions overlap in time, like those indicated in 
Fig. 121 by treating more than one uncontracted photon at a particular time. At present, however, this is beyond reach 
with the computers we have available. On the other hand, reducible multi-photon effects, where the interactions are 
separated in time, as illustrated in Fig. El can be included by repeated use of the procedure described. 

It has been demonstrated that one retarded photon with Coulomb interaction represents by far the dominating 
part of the non-radiative multi-photon exchange for light and medium-heavy elements beyond the standard Coulomb 
correlation (of the order of 99 % for heliumlikc neon) , and the situation can be expected to be similar for the radiative 
part. The small effects due to two irreducible retarded photons without Coulomb interactions can be evaluated 
with standard QED methods, and higher-order effects can with good accuracy be estimated by means of analytical 
approximations. Therefore, it is our belief that the method presented here, when the routines are fully developed, 
should be able to produce accurate results for energy separations, such as the fine-structure separations, for light and 
medium-heavy elements, hopefully down to neutral helium. 
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APPENDIX A: COULOMB-GAUGE INTERACTION 



The terms in the Coulomb-gauge expression (|47|l can be separated into a product of two single-particle potentials 
by the spherical wave expansion, using the relation 



sin(fcri2) 



= kY^{2l + l)ji{kr,)ji{kT2)C\l) ■ C\2), 



1=0 



where C' is a spherical tensor, associated with the spherical harmonics Y^. The Gaunt term then becomes 

oc 

-kY,i2i + i)vl,{kn).vl,{kr,), 

1=0 

where VQ{kri) is 

V^ikr) = aji{kr)CK 

For the scalar-retardation term we use the relation js^ Sect. 5.7], Part II, App. A] 



1 



21 + 1 



I. 



VG + l)(2/ + 3)(f - -)/(r)C:„'+i 
dr r 



d l + l. 



i.i-i 



where C^'^^ is a vector, associated with the vector spherical harmonics 1^^^^, and the relation 

where the left-hand side is a scalar product and the right-hand side a tensor product. Together with [s? 

(.4- - -)ji{kr) = -kji+i{kr) 
dr r 

d ^ [^^^.^^^^^ ^ kji^i{kr) 
dr r 

the final result of the scalar retardation term becomes 



1=0 



where expression for the single-particle potentials for the scalar retardation is written as 

VUkr) = [V{l + mi + 3)ji+i{kr){a&+'y 
+ ^l{2l-l)ji^i{kr){a&-^y . 

The function, /(fc), in the Coulomb gauge expression H47|l then becomes 

1=0 



(Al) 



(A2) 



(A3) 



(A4) 
(A5) 

(A6) 
(A7) 

(A8) 



(A9) 



(AlO) 



APPENDIX B: RULES FOR DIFFERENTIATION 



The difference ratios we use in the formalism presented here are of a special kind and give rise to special handling 
rules (see also ref. 0, App. E). 

If A{£) is an operator function of the (energy) parameter then we define the first-order difference ratio 

SA{£) S££,A{£) Ai£) - A{£') 
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Then 



ob 



The second difference ratio is defined as 

S^A{£) _ Sj,,s„A{£) _ See, Sse^AjS) 



and generally 



It then follows that 



5£ 5£ 



S£ S£ 



Ses" 



S£ 



S£^ 



[Ai£)B{£)] = ^-hl^[A{£)B{£)] = 



Sss' \hs"A{£) 



5£ 
5^A 



5£ 



5£'^ 

B{£") + A{£) 



5£ 5£ 
5££"B{£) 



[A{£)B{£)] 



S£ 



5£'^ 



6£ 



S£ 



S£^ 



(Bl) 



(B2) 



It should be noted that the operator B{£") is unaffected by the differentiation Ssgi. With simplified notations we 
then have 



5'^{AB) _ S^A M m S^B 
~~5£^ ^ ^ ^ j£ j£ ^ J^' 



(B3) 



which can be generalized to 



5£"- ~ ^ j£^ S£("-'"') 



(B4) 



In the case of complete degeneracy we have in first order 



lim 



S££,A{£) dAi£) 



E'^£ S£ 



d£ 



(B5) 



while in second order we have 



Sj,,,„A{£) 1 d'A{£) 



2 d£^ 



(B6) 



and generally 



lim 



5-ee,e„A{£) _ 1 



(B7) 



If we have a product of effective interactions, then they arc separated by model-space states, which generally have 
different energies, 

■ ■ ■ P"V,sP'V,sP, 

where P corresponds to the energy £, P' to £' etc. This appears in multiple folded terms, and the effective interactions 
have the corresponding energy parameter, 

■■■P"V,fi{£')P'V,s{£)P 
If we now form the difference ratio of this product 

±[...p"V,s{£')P'V,i,{£)P], 



implying that the parameter £ is changed, only the last factor is affected, and 

. . P"V,s{£')P'V,^{S)P] ^■■■P"VM£')P' ^-^^j^P 

or generally 

5£ ~ 5£ 
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